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ABSTRACT. Swimming by shape changes at low Reynolds number is widely 
used in biology and understanding how the efficiency of movement depends on 
the geometric pattern of shape changes is important to understand swimming 
of microorganisms and in designing low Reynolds number swimming models. 
The simplest models of shape changes are those that comprise a series of linked 
spheres that can change their separation and/or their size. Herein we compare 
the efficiency of three models in which these modes are used in different ways. 


1. Introduction. Single-cell organisms use a variety of strategies for translocation, 
including crawling, swimming, drifting with the surrounding flow, and others. Some, 
such as bacteria, use flagella, and others, such as paramecia, use cilia to swim, and 
both types use only one mode. However other cells can be more flexible in that 
they either crawl by transient attachments to their surroundings - often called the 
mesenchymal mode, or by shape changes - called the amoeboid mode [13]. The 
former may involve strong adhesion to the substrate or the extracellular matrix 
(ECM) via integrin-mediated adhesion complexes, while the latter depends less on 
force transmission to the ECM or to the surrounding fluid, and instead involves 
shape changes to exploit spaces in the ECM to move through it ( cf. Fig. 1). The 
latter mode can be far more effective and can lead to speeds up to forty times faster 
than those resulting from mesenchymal motion [37]. Cells such as leukocytes, which 
normally use the mesenchymal mode in the ECM, can migrate in vivo in the absence 
of integrins, using a ’flowing and squeezing’ mechanism [25]. While crawling and 
swimming are mechanistically distinct strategies, cells can sense their environment 
and use the most efficient strategy in a given context. 

The spatio-temporal scale of motion of small organisms in viscous fluids fre¬ 
quently leads to low Reynolds number (LRN) flows. Swimmers that use a single 
long, thin flagellum led to the development and application of slender-body theory 
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Figure 1: A sequence of shape changes observed in a Dictyostelium discoideum cell swimming in a fluid 

[9]- 


[27, 20, 22, 16, 11, 23, 21], while microorganisms that swim using a thin layer of 
cilia were first studied by Lighthill for squirming motion of nearly spherical de¬ 
formable bodies [28]. A general review of previous work on swimming appears in 
[26], and here we only analyze models that comprise a number of linked subunits 
and have only a finite number of degrees of freedom. We call these discrete models 
of swimmers. 

Much of the current interest in locomotion at LRN was stimulated by Purcell’s 
description of life at low Reynolds number [35]. In particular, the observation 
that certain classes of shape changes produce no net motion in a viscous fluid 
led to studies on various types of discrete models of swimmers, with the goals of 
understanding how microorganisms swim and facilitating the design of mini-robots 
that swim at LRN. The first discrete LRN model is Purcell's two-hinge swimmer, 
also referred to as Purcell’s three-link swimmer [36] (Fig. 2 (a)). Purcell’s model 
swimmer comprises three connected, rigid segments that are constrained to move in 
a plane and can execute restricted rotations around joints linking the segments. The 
shape is specified by two parameters, the angles between adjacent segments, and 
Purcell showed that one can impose sequences of changes in the angles that produce 
net translation of the swimmer. Despite its geometric simplicity, the relationships 
between geometric parameters, speed and efficiency of swimming are not simple 
[12, 6], but approximations of optimal strokes are known [39]. Various simpler 
linked-sphere models for which both analytical and computational results can be 
obtained have appeared since. The first of these is the Najafi-Golestanian three- 
sphere accordion model (NG) [30, 18, 1] (Fig. 2 (b)), which comprises three rigid 
spheres connected by two slender connecting arms aligned along the ^-direction 
that can stretch and contract in a prescribed form to produce motion. Since the 
forces that expand or contract the arms are directed along them it can only result 
in translation - it never rotates. Another linked-sphere model is the pushmepullyou 
swimmer (PMPY) [7] (Fig. 2 (c)), in which two spheres that can expand or contract 
radially are connected by an extensible arm. Analytical and numerical studies of 
the NG and PMPY models have been done heretofore, and their efficiency and 
the optimality of various strokes have been investigated [2, 3, 4]. Recently we 
have analyzed a three-sphere volume-exchange or breather model (VE) in which 
the spheres are linked by rigid connectors but exchange volume [42] (Fig. 2 (d)), 
the details of which will be discussed in Section 3. 




CELL SWIMMING 


3 


(a) Purcell 


\ 


/ 

/ 

V 

\ 


/ 

/ 

\ 

\ 

/ 


(b) NG 


(c) PMPY 



(d) VE 




Figure 2: Low Reynolds number swimming models: (a) Purcell's 3-link swimmer [36]; (b) Najafi- 
Golestanian’s 3-sphere model (NG) [31]; (c) Pushmepullyou (PMPY) [7]; (d) The 3-sphere volume- 
exchange model (VE) [42], 


A central problem in the analysis of both biological LRN swimmers and mini¬ 
robots is whether a cyclic sequence of deformations results in significant movement, 
and if such a swimming mode is efficient by some measure. This has been studied for 
Purcell’s swimmer [39], the PMPY swimmer [5] and cilia-based swimming [29, 32], 
Lighthill’s definition [28] provides one metric of efficiency and several others have 
been used, but we will introduce a new criterion to measure the performance of 
LRN swimmers. 

Of course in reality the various shape changes that have been analyzed require 
internal forces that generate the shape changes needed for propulsion of the or¬ 
ganism. In biological organisms this interior problem involves the biochemical and 
biophysical changes in the cytoskeleton needed to produce the necessary intracellu¬ 
lar forces and shape changes, but here we simply prescribe the shape changes and 
treat the exterior problem. An integrated model that includes sensing the environ¬ 
ment and controlling the shape changes so as to move is still beyond reach for even 
a single-cell organism. 

2. Movement by shape changes — the exterior problem. The Navier-Stokes 
equations for an incompressible fluid of density p, viscosity p, and velocity u are 
du 

P l)i + ^ = V ’ a + f ext = -Vp + pA U + f ext , (1) 

V • u = 0 (2) 

where er = -p^ + /i(Vu+(Vtt) :r ) is the Cauchy stress tensor and f ex t is the external 
force field. Herein we assume that the swimmer is self-propelled and does not rely 
on any exterior force, and therefore we require that f ex t = 0. The Reynolds number 
based on a characteristic length scale L and speed scale U is Re = pLU /p, and 
when converted to dimensionless form and the symbols re-defined, the equations 
read 

du 

ReSl— —I -Re(u-X7)u = —Vp + Ait, 
at 


V • u 


0 . 


( 3 ) 
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Here SI = uiL/U is the Strouhal number and ui is a characteristic frequency of 
the shape changes. When Re<C 1 the convective momentum term in (3) can be 
neglected, but the time variation requires that ReSl = uL 2 /v -C 1. When both 
terms are neglected, which we assume throughout, the flow is governed by the Stokes 
equations 

/tA u — Vp = 0, V • u = 0. (4) 

We also only consider the propulsion problem in an infinite domain and impose the 
condition ulx-^oo = 0 on the velocity. 

In the LRN regime time does not appear explicitly, momentum is assumed to 
equilibrate instantaneously, and bodies move by exploiting the viscous resistance 
of the fluid. As a result, time-reversible deformations produce no motion, which is 
the content of the ‘scallop theorem’ [35]. In the absence of external forces due to 
boundaries or other fields there is no net force or torque on a self-propelled swimmer 
in the Stokes regime, and therefore movement is a purely geometric process: the net 
displacement of a swimmer during a stroke is independent of the rate at which the 
stroke is executed, as long as the Reynolds number remains small enough. Amoebae 
of the slime mold Dictyostelium discoideum have a typical length L ~ 25/rm and can 
swim at U ~ 3/xm/min [40]. Assuming the medium is water (p ~ 10 3 kg m , // ~ 
10~ 3 Pa • s), and the deformation frequency u ~ 1/s, Re ~ 0(1O -6 ) and SI ~ 
0(1O -4 ). In fact the experiments are done in oil that is significantly more viscous 
[9], and for similar cells one can neglect both inertial terms. 

Suppose that a swimmer occupies the closed compact domain f l(t) C R n (n = 
2,3), at time t, and let dfl(t) denote its prescribed time-dependent boundary. A 
swimming stroke T is specified by a time-dependent sequence of the boundary l(t), 
and it is cyclic if the initial and final shapes are identical, i.e., c>f2(0) = dli{T) 
where T is the period [38]. The swimmers’ boundary velocity V relative to fixed 
coordinates can be written as a part v that defines the intrinsic shape deformations, 
and a rigid motion U. If u denotes the velocity field in the fluid exterior to fl, then a 
standard LRN self-propulsion problem is : given a cyclic shape deformation specified 
by v, solve the Stokes equations subject to 

/ er ■ n = 0, / r A (er ■ n) = 0, u| on(t ) = V = v + U, u| x -^oo = 0 (5) 

JdQ(t) JdQ.(t) 

where n is the exterior normal, and the integrals are the force- and torque-free 
conditions. 

In order to treat general shape changes of a cell defined by fl(t.) £ R 3 with 
boundary <9f1(f), one must solve the exterior Stokes equations (4) for it, with a 
prescribed velocity v(t) on <9f l(t) and subject to the decay conditions u ~ 1/r and 
p ~ 1/r 2 for r —^ oo. It is known that the solution has the representation 

«(*) = -— 1 — [ G(x, y) ■ f(y)dS(y) + / v(y) ■ T(y, x) ■ n dS{y) (6) 

Xirp Jon(t) o' 71 " JdQ(t) 

where G is the free-space Green’s function, T is the associated third-rank stress 
tensor, and f = er ■ n is the force on the boundary [34]. The constraints that the 
total force and the total torque vanish determine the center-of-mass translational 
and angular velocities. When x £ dQ(t) this is an integral equation for the force 
distribution on the boundary, the solution of which determines the forces needed 
to produce the prescribed shape changes. The approach is similar in spirit to what 
has been done for cells crawling on a deformable substrate, where the substrate 
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deformations are given and the forces exerted by the cell are the solution of a 
Fredholm integral equation [8, 14]. 

The free space Green’s function or Stokeslet has the form 


G(x,x 0 ) 



( 7 ) 


where I is the unit second-rank tensor, r = x — xo, and r 
velocity field generated by a point force f at the origin is 


u{x) 


G(x, 0) f 

8nfi 


|x — xq|. Thus the 


( 8 ) 


G(x, 0)/(87iyz) is called the Oseen tensor. Three other basic solutions that are 
needed are those for a rigid sphere pulled through a quiescent fluid, for a radially 
expanding or contracting sphere, and for the interaction between two spheres. 

When a sphere of radius a is pulled through a quiescent fluid with a steady force 
F under no-slip conditions at the surface, the resulting flow field is given by 


u(r) = F ■ (1 + 


« 2 V 2 . G(x,a: s ) 

6 87 TfJ 


F 

87r/ir 


, rr a 2 

r T rri 

I T —pr T —— 

I — 3 —^ 

r 2 3r 2 

r 2 J 


( 9 ) 


where x s is position of the center of the sphere and r = x — x s . The second 
term represents the degenerate quadrupole needed to satisfy the no-slip boundary 
condition at r = a, but it is small when a/r <C 1. The resulting velocity of the 
sphere is given by Stoke’s law [24] 


F = 6irnaU. 


( 10 ) 


This can be obtained directly from (6) by setting v = 0 in the second integral, 
expanding the Green’s function, and defining the total force on the sphere as 


F = [ {(y)dS(y). 
Jen 


( 11 ) 


A second basic solution is the velocity field u produced by a radially expanding 
sphere, which can be generated by a point source at the center x s of the sphere [34]. 
The corresponding velocity is 

u = a-^ (12) 

where r = x — x s and a is a constant to be determined. The no-slip boundary 
condition at the surface r = a requires that 


u(r) 



da r 
d t a 


and therefore a = da 2 , and 



( 13 ) 


(14) 


where v = 47ra 3 /3 is the volume of the sphere and r = r/r. By combining equa¬ 
tions (9, 13) we obtain the velocity field for the combination of the pulled and 
expanding sphere, namely [7] 


u (r; a, F, d) = 


1 


247 Tfj,r . 


(3+-j)F + 3(l--2)(F.r)r 


47r r 2 


r. 


(15) 
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The last basic solution needed involves the interaction between two spheres 1 . 
Suppose that the *th (i = 1,2) sphere has radius a,i(t), is centered at x,(t) and is 
subjected to a drag force F i{t) due to its motion. The translational velocity of the 
itli sphere consists of two parts: Up 0 = (67r/xai) _1 Fj that results from the drag 
force Fj exerted on i, and the other a perturbation part <5U,; that is due to the flow 
generated by the other sphere. In particular, since the other sphere is translating 
and expanding, i)U,; can be further decomposed into two parts: <5U* due to the 
translation of the other sphere, and <5U® that results from its radial expansion. 
Hence we have the following decomposition of the translational velocity of the zth 
sphere. 


U,= 


F; 


6nfMi 


- <5U- + <5U? 


(16) 


Here <5U* arises from a flow given by equation (9) and is given by 


SVl = 


(l + |v 2 )u(r) 


1 


87 Tfll . 


1 + 


r=x i -x i 

af + a 2 


3l 2 


L M*- 


a 2 \ 11 


l 2 )l 2 


(17) 


where 1 = x, — Xj and l = |x, — x, [10]. The velocity S Uf is resulted from a flow 
given by equation (13): 

,2 


SUf = 


(l + ~g~V 2 ) u (r) 


r—Xi —x,- 


= u( Xi - Xj) = 


l 3 


(18) 


Altogether, <5U- and <SUf induce a higher-order perturbation in U, but as we shall 
see, these terms are neglected in the existing analyses of linked-spheres. 

Next we use these solutions in the analysis of various models, and we begin with 
the pure volume-exchange (VE) model. 


3. The 3-sphere volume-exchange model. Some cells produce membrane pro¬ 
trusions called blebs that emerge when the membrane detaches from the cortex 
locally and the excess internal pressure forces fluid into the bleb [17, 33]. When this 
occurs repeatedly over a cell surface, as in Fig. 3(a), it may result in an oscillatory 
motion of the cell. In Fig. 3(a) the cell blebs blebs profusely with little net transla¬ 
tion, whereas Fig. 3(b) shows a motile, blebbing Dictyostelium discoideum cell. If 
bleb formation is restricted to the leading edge as in (b), forward motion is driven 
by contraction of the cortical network at the rear of the cell. In either case one 
can understand the dynamics in terms of mass or volume exchange between differ¬ 
ent parts of the cell. The protrusions are usually approximately hemispherical and 
thus a linked-sphere model may be a good choice for a study of blebbing. However 
most existing linked-spheres models require significant changes in the length of the 
connecting links, which is not realistic in blebbing cells. This led us to develop a 
model that better describes blebbing dynamics [42]. 

In a minimal VE model there are three spheres linked by two rigid, mass-less arms 
of fixed length l (Fig. 2(d)). A sphere can only expand or contract in the radial 
direction (i.e., = Gq(i)), and can only exchange mass with its neighbor(s), i.e., 1 

with 2 and 2 with 3 but not 1 with 3. There is little evidence that swimming cells 
exchange significant material with the surrounding fluid, and thus we impose mass 


1 For simplicity we assume that the interactions in a general configuration of spheres are pairwise 

additive. 
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Figure 3: (a) Blebbing on a melanoma cell: myosin (green) localizes under the blebbing membrane 

(red) (from[15]) (b) The actin cortex of a blebbing Dd cell migrating to the lower right. Arrowheads 
indicate the successive blebs and arcs of the actin cortex (from[43]). 


conservation on the ensemble. We are primarily concerned with the fluid-structure 
interaction and for simplicity we assume that the density in all three spheres is 
the same constant, in which case the mass conservation is equivalent to volume 
conservation. We also ignore dissipation within cells, and with these constraints, it 
is easily seen from the “scallop theorem” [35] that a cyclic two-sphere model cannot 
swim, since it has only one degree of freedom. Therefore a minimal model must 
comprise at least three spheres. 

The velocity field is given by equation (15), and to obtain asymptotic solutions, 
we assume that the distance l between either pair of neighboring spheres is larger 
than the radii cq (i = 1,2,3), i.e., ai/l <C 1, which simplifies the computation yet 
captures basic aspects of the movement of blebbing cells. The VE model is linear, 
thus the velocity U, of the ith sphere and the force F, exerted on the ith sphere 
are all along the x-direction. 

Ui = UiQ x , F i = Fig x 

Because the Stokes problem is linear, the velocity Ui of each sphere is given by 
equations (16, 17, 18), wherein we only retain the leading order terms in the per¬ 
turbations <5U* and <5U f. The asymptotic solution for the Ui is given by 


Ui r 

Fi 


/ f 2 

V2 \ 

+ ( 

■ f 3 

V2 \ 

(19) 

Gn^iai 


V 47t fil 

47 tP ) 


16irl' 2 ) 

U 2 r 

f 2 

+ 1 

( Fl 

1 7)1 "l 

+ ( 

■ f 3 

i’3 \ 

( 20 ) 

677 fia 2 

V47 Tfll 

4tt( 2 / 

.47 T/xl 

47 tI 2 ) 

U 3 r 

f 3 


f ^ 

1 ^ ' 

h 

/ f 2 

1 i2 "l 

( 21 ) 

6n^a 3 


K8n[il 

167r ( 2 

\47 Tfll 

47 tP ) 


where Vi = 47raf/3 and 1)2 = —v 1 — i> 3 . Since the connecting arms have fixed length 
l, U\ = U 2 = U 3 , which defines the translational velocity U of the swimmer. Since 
the swimmer is linear it is necessarily torque-free, and the force-free constraint is 


F\ + F 2 + F 3 = 0 . 


(22) 
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The volume conservation constraint reads 


iii + i>2 + i> 3 = 0 


(23) 


and (19 — 23) lead to the following asymptotic solution for the swimming velocity 
of the model [42], 


u _ (ai+a 2 - |a 3 )t>i - (a 3 + a 2 - f ai)t > 3 
47 r/ 2 (ai + a 2 + a 3 ) 


(24) 


Next, we consider the power P = fj" f (t)U(t)dt required to propel the swimmer. 
The stress on the surface of the expanding sphere is a = pv/(na 3 ) [7], and therefore 
the power required to expand one sphere is 


. 9 • 4 /r - 2 

47 ra aa = av = —v 
6v 


(25) 


Therefore the total instantaneous power expended by the swimmer in transferring 
volumes between the spheres is 


P = 


4/r 


p? , v 2 , * 3 2 1 

_ M 

\( 1 , 1 

. 2 2 . . 

( 1 , 1 \. 2 l 

— + — + — 


( ^3 H 3 

)v l + -^ViV 3 + 

“3 H 3 ) ^3 

L'Ci v 2 n 3 J 

7 r 

VaJ dj 

a 2 

\(l2 «3 / 


(26) 


Finally we define the performance of a stroke T as the ratio of the translation 
per cycle to the energy expended in a cycle. 


I Ip U(t)dt\ 

fo P(t)dt 


(27) 


This has units of 1 /force. 

When the volume changes are small several conclusions can be reached analyt¬ 
ically [42, 41]. Since a swimming stroke is a closed path in the V\ — v 3 plane or 
equivalently, a closed path in the a 3 — a 3 plane, we find the following relation be¬ 
tween the differential displacement da; and the differential controls (dai.da 3 ) from 
(24). 


~ dX = ^ 




03 


7 


4 ai + a 2 + a 3 


)dai — a§(l— - 


a 1 


4 a\ + a 2 + a 3 


)da 3 


(28) 


Here dX > 0 represents an infinitesimal displacement to the right in Fig. 2(d). The 
bar in dX indicates that the differential displacement is not an exact differential. 

To determine the direction of swimming, note that from equation (28) we may 
assume, without loss of generality, that da 3 = 0 and dai > 0 , which means that 
sphere 3 does not change, and sphere 1 is expanding while sphere 2 is contracting. 
For a 2 large enough so that 1 — 7 a 3 /(ai+a 2 +a 3 ) > 0 always holds, we have dA' > 0, 
which means that the swimming direction is from sphere 1 to sphere 2. Hence we 
have the following conclusion, which also applies to the PMPY swimmer. 


Conclusion 1. When only one pair of adjacent spheres is involved in volume ex¬ 
change, and when the central sphere is large enough, the direction of swimming is 
always from the expanding sphere to the contracting one. 


Next, using Stokes’ theorem, the translation SX corresponding to an infinitesimal 
closed loop is 


4 1 2 


ajda 


a 3 


a.\ T a 2 Y a 3 


and. 


a 1 


3^01 


a l + a 2 + a 3 


daiA da 3 


(29) 
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where dai A da 3 denotes the signed area enclosed by the loop. From this one can 
show the following. 


Conclusion 2. For strokes such that F is homotopic to the unit circle, increasing 
the stroke amplitude will increase the net translation of the stroke, while increasing 
the initial radius 020 of the central sphere (with aio and a 3 o unchanged) will decrease 
the net translation. Moreover, we have the approximation 

|X(r)| ~ jArea(n) 

where e ~ ai/l and Q is the region enclosed by T and Area(f2) is its signed area. 
The proof of this is given in [41]. 


3.1. Numerical computations. Next we prescribe cycles of shape changes in the 
controls (di,a 3 ) and compute the displacement and performance measure numeri¬ 
cally. In particular, we investigate how the following characteristics of the system 
affect the net translation X = f Q U(t)dt and the performance e of the swimmer 
after a full cycle (T = 1). 

1. L, which measures the distance between a pair of neighboring sphere; 

2 . tt, r 2 , which measure the amplitude of shape deformations; 

3. s, which measures the size of the central sphere. 

We use the following protocol in varying these parameters. 


Arm length (fixed) : 

Controls: 

Radius of the central sphere: 
We consider three cases. 


l\ — I2 — L 

ai(f) = Ro +r\ cos 27ri, 03 (f) = Ro + sin27rt, 

a 2 ( 0 ) = s, a 2 (i) = (Ay to t - ai(t ) 3 ^ o 3 (f ) 3 ) 3 


I First, we set n = r 3 = 1, -Ro = 2 and s = 3, and test different values of L, i.e., 
arm length (Fig. 4(a,b)). Fig. 4(c) gives the initial profile of the swimmers with 
L = 10 and L = 50, respectively. Fig. 4(a,b) show that both the translation X 
and the performance e decrease as the arm length increases. That is to say, with a 
fixed amount of body mass, a too long body is not a good strategy for swimming. 
The effect on the translation X can be seen from equation (24) and Conclusion 2, 
with a fixed stroke rending the same Area(f2), the translation X scales in the order 
of eF 1 II = aF 2 . Hence as the arm length L increases, A' decreases quickly. As for 
performance e (equation (27)), it is the ratio of translation to energy over a cycle, 
with power P given by equation (26), from which we clearly see that the arm length 
L does not enter into the expression of P. Hence e decays similarly to A', namely, in 
the order of l~ 2 . 


II Next, we set L = 50, s = 10, Ro = 10, and test different values of ri and r 3 
(Fig. 4(c,d)), i.e., the stroke amplitude. Fig. 4(e) gives the initial profile of the 
swimmers with r 3 = 1 always, but ri = 1 or n = 9, respectively. From Fig. 4(c) 
we see that The translation increases as either n increases, but if one of ri is small, 
the increase of translation due to the other r, is small. This can be explained by 
Conclusion 2, in general a large ri or/and r 3 indicate a large stroke, i.e., a large 
Area(ff), which clearly induces a large translation X — except that when one of ri is 
small, then no matter how large the other ri is, we have Area(f2) ~ 0, which results 
in a A ~ 0. As for performance, it is difficult to obtain a general analysis between 
e and Area(fi), yet Fig. 4(c) indicates that e has the similar behavior as X (though 
maybe not to the same order), i.e., e increases as either n increases, if the other n 
is not too small. 
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Figure 4: (a,b) Translation X and performance e as a function of the fixed arm length L. (c) The 
initial shape of the swimmer when L = 10 or L = 50 in simulations shown in a& b. (d,e) Translation X 
and performance e as a function of the stroke amplitude (t*i,t* 3 ). (f) The initial shape of the swimmer 
when 7*1 = 1 or r± = 9 in simulations shown in (d,e). Notice that although both r\ and 7*3 determine 
the stroke amplitude, the initial shape is only related to r±. (g,h) Translation X and performance e in 
relation of the (initial) size of the central sphere s. (i) The initial shape of the swimmer when s = 3 
or s = 15 in simulations shown in (g,h). The scales are the same for panels (c,f,i). 


Ill Finally we set L = 25, n = r 3 = 1 and Ro = 2, and test different values of s 
(Fig. 4(e,f)), i.e., the (relative) size of the central sphere to the side spheres. Here 
we would like to point out such a test can be considered a complement to the test of 
stroke amplitude. Consider that we may scale the total volume of the whole object 
to be the same with respect to different values of s, then a model with a small central 
sphere (Fig. 4(i), s = 3) is translated to one with a large stroke, and on the other 
hand a large central sphere (Fig. 4(i), s = 15) corresponds to a swimmer with a small 
stroke. Hence the results can be predicted from previous discussion and are shown in 
Fig. 4(e,f): both the translation and the performance decrease as the initial radius 
of the central sphere is increased. Yet we illustrate such behavior in regard to s. 
Observe equation (29) and we find that with both the stroke amplitude dai A da 3 
and the arm length L fixed, X is proportional to the following quantity: 


ajda 


03 


ai + (22 + &3 


+ Oj^3a\ 


CLl 


&l ~\~ CL2 G-3 


(30) 


which clearly decreases as the central sphere gets bigger. On the other hand, equa¬ 
tion (26) indicates that a bigger central sphere results in small power P. To evaluate 
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the effect of the size of the central sphere on the performance e, we conduct an as¬ 
ymptotic analysis with the scenario 01,03 <C 02 . In such a case, the quantity given 
by equation (30) approximates 0, then equation (29) shows that 5X ~ 0. However 
from equation (26) we obtain the following behavior P when ai, 03 -C 02 : 


M T 1 -2 , 1 -2l 

_ ~3 V 1 H- 3 V 3 

7T 1 O^ dj J 


which in general does not vanish. Hence as the ratio of average translation to average 
power, the performance e vanishes in the scenario 01,03 <C 02 . 


In conclusion, in order to reach longer net translation X or better performance 
e, the 3-spliere volume-exchange model should be designed so that the connecting 
arm is short, more mass are exchanged among the spheres though the cycle, and the 
central sphere should not be big comparing to the two side spheres. 


4. A comparison of the three linked-sphere swimmers. Here we compare 
the performance of the three linked-sphere swimmers: the NG 3-sphere model, the 
pushmepullyou (PMPY) and the VE model. A summary of the analysis that leads 
to the velocity and power of the NG and PMPY models can be found in Appendix 
A and in [31, 19, 7]. To standardize the comparison between them, we stipulate 
that the total volume in all spheres are the same for each model, and the stroke 
amplitudes are the same. We prescribe strokes for each model as follows. 

• NG: 

Ri = R2 = R3 = 2 rG 

h(t) = L + cos(27r t), I 2 {t) = L + sin(27rf) 

• Pushmepullyou: 

Ri = 2rp + cos(2-7t t) 

i?. 2 (0) = 2r P , R\ + Rl = ^-V tot 

4-7T 

l(t) = L + sin(2-7rf) 

• Volume-exchange (VE): 

li — I 2 = L 

Ri(t) = 2 rv + cos(27rf), -fib(t) = 2 ? - vb + sin(27rf) 
i? 2 (0)=2r v , Rl + Rl + Rl = i-Hot 

47T 

where scales rQ,rp,ry are chosen so that the total volume of all spheres in each 
swimmer are the same. Without loss of generality, we may choose 7'q = 1. 

First, the velocity U(t) and the power P(t) within a cycle that result from the 
above prescribed strokes are given in Fig.s 5(a — d), for L = 6 or 30, respectively. 

For translation, comparing Fig. 5(a) and (c) we see that (7(f) for PMPY does not 
change much, while it almost vanishes for NG and VE when L = 30. On the other 
hand, P(t) in Fig. 5(b) and Fig. 5(d) are quite much similar. 

Moreover, we find that the PMPY model is clearly superior to the others (Fig. 5(e,f )), 
in regard to both translation and performance over the entire range of L, even when 
L —> 00 . The reason can be found in the asymptotic solutions of U and P for the 
three swimmers. While the leading term of U for PMPY is 0(1), it is C7(f -1 ) and 
0(l~ 2 ) for NG and VE, respectively. The leading order term in the power P is 
0(1) (equation (26)) for all three swimmers, hence the arm length L does not have 
much influence on the power or the performance e. The leading order term of e is 
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(a) Velocity U(t) (b) Power P(t) 



(c) Velocity U(t) (d) Power P(t) 



L L 


Figure 5: A comparison of the three linked-sphere swimmers. (a,b) U(t) and P(t) of the swimmers 
within a cycle with L = 6. (c,d) U(t) and P(t) of the swimmers within a cycle with L = 30. (e,f) Net 
translation X and performance e of the three swimmers with different values of L. 


the same as that of U, i.e., 0(1), 0(1 x ),0(l 2 ) for PMPY, Golestanian and VE, 
respectively. 

5. Mixed controls result in more net translation and better performance. 

All three models have two degrees of freedom, which are of two types - a change in 
the arm length l, or a change in the sphere radius a. Different combinations of the 
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controls result in different swimming behaviors, and from the results above we find 
that the best choice is the mixed strategy, i.e., (I, R), which is adopted by PMPY. 
A combination of the same kind - (h,h) or (i?i,i? 2 ) — is not advantageous, and 
they produce comparable net translation and performance (Fig. 5). 

The priority of the mixed strategy is based on the following two principles: 

1 . I results in a velocity with leading order 0 ( 1 ), while the leading order that 
results from R is only 0(l~ 2 ). Hence to increase the net translation, one 
should incorporate a change in arm length. 

2. The leading order terms in the velocity U(t) should not be, or even approxi¬ 
mate, an exact differential, because they then make no or little contribution 
to the net translation. 


We have discussed the first principle above, which clearly explains why PMPY 
is better than the other two swimmers. However this does not explain why the net 
translation behaves similar for NG and VE, with leading order terms that scale like 
e>(Z -1 ) and 0(l~ 2 ), respectively. In fact, it is only in the case that all three spheres 
are of equal size that the velocity scales like 0(/ _1 ) for the NG swimmer. When 
tspheres are of different sizes are involved, U(t) will actually scale as 0(1), i.e., the 
same as U(t) in the PMPY model. However, the net translation X still turns out 
to be of order 0(l ~~ 2 ), which is the same as for the VE swimmer. 

To understand this, we analyze the asymptotic solution of U(t) for an NG swim¬ 
mer which is approximated to the order 0(/ _1 ) (equation (Al), [19]). The equation 
is complex, but it can be written in the following form. 


U — A(di + B 0 I 2 + 


4 ^ , 4 

h 


(2) 


+ 


A 


( 12 ) 


h + h 


)h+( 


B 


(i) 


h 


+ 


B 


( 2 ) 




( 12 ) 


h + h 


h+0(-)( 31) 


where all coefficients A 4 , B^" 1 are functions of ai, 02,113 only and do not depend 
on li, I 2 or time t. 

The leading order term of U, denoted as f7( 0 ), is the combination 


U(o) — AqIi + B 0 I 2 


and the integral over a whole cycle gives 


*( 0 ) = 


r T r T . . 

/ t/( 0 ) d t= (Aoh + B 0 I 2 ) dt = Ao/i | t _ 0 +-Bo4 | t=0 =0. (32) 

Jo Jo 


Next, the 0(1 x ) term of U, which we denote t/(i), is given by the following. 


U( i) = 


/A (1) 

( — 
V h 


A\ 


( 2 ) 


A 


( 12 ) 


h + h 


h + 


B 

\ h l 


( 2 ) 


+ 


B 


( 12 ) 


h + h 


In general, the integral f 0 U\ df does not vanish, but we have the relation 

h(t) = L + Sk(t) 


where L is the fixed part and 5U is the deformation part. When L is sufficiently large, 
so as to ensure that the higher-order interactions between spheres are negligible, we 
have that Sli <C L , and thus U(i) can be approximated as 


U(i) 





( 12 ) 


2 L 


)h+( 


B (J) ^ b (2) 


B \ 


( 12 ) 


2 L 
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Thus again, 

A'(i) 


From equations (32,33), we see that although in general the 0(1) and 0(l~ x ) 
terms do not vanish in U(t), the 0(1) term is an exact differential and the 0(1 -1 ) 
term is approximately an exact differential, and hence the net translation X ~ 
0(l~ 2 ). This is the same as the leading order that results from varying the radii in 
the VE model. 

Fig. 6(a,b) give U(t) within one cycle for four NG swimmers, whose radii are given 
in Fig. 6(c), with L = 10 or L = 100. Again, to make a fair comparison we require 
the total volume of all three spheres in each swimmer are the same. We do not 
use equation (Al) in [19] to solve for U(t), instead we numerically solve the whole 
system (equations (34 - 36)). From Fig. 6(a,b) we see that for the same swimmer, 
the amplitude of U(t) within a cycle is almost of the same scale when L = 10 or 
100, yet the net translation is very small with either value of L. However we do 
observe that among different choices of the sphere sizes, net translation favors the 
equal sized spheres (SO) - i.e., SO results in the most net translation yet it requires 
the least amplitude of U(t ) among the four swimmers. 

On the other hand, for PMPY the leading term in U ( t ) is 


U{ i) d t 


4 1} A[ 


( 2 ) 


,( 12 ) 


L 


+ 


2 L 


\ f /. r\i _L 1 

/ 

b[ 2) 

V 2 ’ 

) [ 

)Jo 11 

l L 

' L 

2 L 



h df 


(33) 


U( o) = 


ui — a 2 
2 (cii + 02) 


i ~ o(i) 


which is not an exact differential, and as a result, the leading order term in the net 
translation gives 


A'( 0 ) = 


Cll — a 2 


'0 2 (ai + <22) 


i At ~ 0(1) 


in general, which explains the better performance of the PMPY swimmer as com¬ 
pared with the Golestanian or VE swimmer. 


6. Discussion. We have compared the three most widely-studied discrete swim¬ 
mers and have shown that the PMPY swimmer performs best under the imposed 
conditions. This conclusion is important for the design of mini-robots, but of course 
real systems are more complex, and as we indicated earlier, the VE model is a more 
realistic model of cellular motion. An analysis of continuum models of swimmers 
that more accurately reflect actual shape changes will be reported elsewhere. 


Appendix A. A summary of the NG and PMPY models. The NG swimmer 
(Fig. 2(b)) consists of three spheres with radii a* (i = 1,2,3) and two connecting 
arms with length li(t) (i = 1,2) [30, 18, 1]. While the spheres are rigid, the con¬ 
necting arms can stretch or contract. In the case when ai/l <C 1, the velocities of 
the spheres (Ui) are related to the forces exerted on the spheres (Fj) via the Oseen 
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(c) Radii of sphere in each swimmer 



R1 

R2 

R3 

SO: 

3.00 

3.00 

3.00 

SI: 

3.40 

1.34 

3.40 

S2: 

2.00 

4.02 

2.00 

S3: 

2.00 

2.08 

4.00 


(d) Cartoon illustration of each swimmer 


SO 

o 

SI 

o 

S2 

+- 

S3 

o 






Figure 6: (a) U(t) of four Golestanian swimmers, with L = 10. (b) U(t) of four Golestanian swimmers, 
whose spheres have correspondingly the same size as in (a),but with L = 10. (c) The radius of each 
sphere in each swimmer of (a,b). (d) Initial shapes of the four swimmer (with L = 15 for better 
visualization effects). 


tensor: 


The velocities 


The system is 


and equations 
defined as 


Ux 

F\ F 2 F 3 

= + + 

6n^iai Anfxli 47t/i(/i + l 2 ) 

(34) 

u 2 

-Fi F 2 F 3 

Anfili 6nfj,a 2 ^tt^1 2 

(35) 

U 3 

_ Fi F 2 F 3 

47 + l 2 ) 4tt/i 1 2 Qitfm 3 

(36) 

are related 

via the following relations: 


u 2 

II 

1 

II 

(37) 

force-free 

F\ + F 2 + F 3 = 0 

(38) 


(34 - 38) are a closed system. The velocity of the whole swimmer is 


u = |(C/i + u 2 + u 3 ) 


( 39 ) 
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In the case that the spheres are of equal size, we have the asymptotic solution 


U = - 


a [l 2 — li 




(40) 


6 LZi +12 

The power consumption of the swimmer comes from dragging the spheres, thus 

P = JTjt/i + F 2 U 2 + P 3 P 3 (41) 

and in the case that the spheres are of equal size the above equation simplifies to 


47 T/m 


a a 

l\ 2l 2 


h + h 


ii + 

5 a 


' a a 

2li 2l 2 2(l\ + 1 2 ) 


a a 

2l\ l 2 

i\l 2 . 


h + h 


(42) 

(43) 


The PMPY swimmer (Fig. 2(c)) consists of two spheres with radii dj(t) (i = 1, 2) 
and one connecting arms with length l(t). The spheres can expand or contract 
in the radial direction, and the connecting arms can stretch or contract. When 
ai/l < 1 , the velocities of the spheres ([/*) are related to the forces exerted on the 
spheres (Pi) via the Oseen tensor 


E/i = 

U 2 = 


Fi 

6TTfMLl 

F 2 


l 2 


ai 


:a 2 . 


671/102 l 2 

The velocities are related via the following relation. 

U 2 - I/-, = i 

Again, the system is force-free 

Pi + P2 = 0 

and the total volume of the two spheres is conserved 

a^di + a\a2 = 0. 

Equations (44 - 48) are a closed system and the velocity of the swimmer is 


„ 1 , TT TT . ai - a 2 ■ Oi . 

u = r(E/i + u 2 ) = l + 

2 2(ai + CI 2 ) 


1 2 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


The power consumption P(t) of the swimmer comprises two parts: Pdrag that results 
from the drag force on the spheres, which is given by 


Fdrag — F1U1 + F 2 U 2 

and P exp that results from the radial expansion of the swimmers 

P ex p = 167 T/l(aidi + 0202). 

Hence the power expended is given by 

/I 1 \ -i. 

P = Qnfi( - 1 -) l 2 + IQt: + a 2 a\). 

\di &2' 


(50) 

(51) 

(52) 
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